Giant Anisotropic Magneto-Resistance in ferromagnetic atomic contacts 
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Magneto-resistance is a physical effect of great fundamental and industrial interest since it is 
the basis for the magnetic field sensors used in computer read-heads and Magnetic Random Access 
Memories. As device dimensions are reduced, some important physical length scales for magnetism 
and electrical transport will soon be attained. Ultimately, there is a strong need to know if the 
physical phenomena responsible for magneto-resistance still hold at the atomic scale. Here, we 
show that the anisotropy of magneto-resistance is greatly enhanced in atomic size constrictions. We 
explain this physical effect by a change in the electronic density of states in the junction when the 
magnetisation is rotated, as supported by our ab-initio calculations. This stems from the "spin- 
orbit coupling" mechanism linking the shape of the orbitals with the spin direction. This sensitively 
affects the conductance of atomic contacts which is determined by the overlap of the valence orbitals. 

PACS numbers: PACS numbers: 75.70.Kw, 72.75.Gd 



The effect of an external field on the resistivity of pure 
ferromagnetic metals (the magnetoresistance-MR) was 
the subject of intense research work in the second half of 
the 20ieth century. The field has seen a renewed interest 
in the past fifteen years with the discovery of giant effects 
in systems combining magnetic and non-magnetic mate- 
rials. This Giant Magneto-Resistance (GMR) has had a 
tremendous impact both through its industrial applica- 
tions as read-heads and Magnetic Random Access Mem- 
ories as well as for triggering the field of " spintronics" 
[1] , aiming to use the spin of the charge carriers in elec- 
tronic devices with higher functionalities. As the pressure 
towards miniaturization increases, it is important to un- 
derstand how magnetoresistive effects are influenced by 
size reduction. In constrictions of dimensions close to 
the Fermi wavelength, boundary conditions enforce that 
transverse electronic modes are quantized which results 
in the discreteness of propagating electron modes. 2-D 
electron gases are archetypical systems in which the con- 
ductance is quantized in units of 2e^//i. In metals where 
the Fermi wavelength is typically 3 A, one needs to reach 
atomic dimensions in order to observe such effects [2]. 
But because even in the single atomic regime several or- 
bitals overlap, one normally finds that several conduc- 
tion channels are opened with imperfect transparency, 
i.e. each channel has a transmittance associated to it (a 
coefficient between and 1). Calculations seem to indi- 
cate that 4 or 5 channels participate significantly to the 
conduction of 3d transition metal nanocontacts [3, 4]. 
The magneto-resistance obtained when one side of the 
contact flips its magnetization is enhanced compared to 
that in the bulk [4, 5] and values of the order of 20% 
have been reported in some careful experiments [6, 7]. 
On the other hand, one could expect some dependence of 
the conductance to the direction of the magnetization be- 
cause changing the spins' direction will affect the orbitals 
through a mechanism known as spin-orbit coupling. One 



can describe this interaction by a term in the system's 
Hamiltonian written XL.S where A is called the spin-orbit 
constant, L and S are the orbital momentum and spin 
vectors. In a solid, molecular orbitals lose most of their 
angular momentum, an effect known as " quenching" . As 
a result, the resistance variation with the angle between 
the local moments and the electrical current lines: the 
Anisotropic Magneto Resistance (AMR), is only at most 
a few percent. This is actually the oldest known magnetic 
effect on electronic transport in ferromagnets, which was 
discovered in 1857 by W. Thomson [8]. When the dimen- 
sionality of the system is reduced, like in mono-atomic 
wires, orbital moments are much larger than in the bulk 
[9]. One can then expect spin-orbit mechanisms, like the 
AMR, to be enhanced. Unfortunately, it is extremely dif- 
ficult to study electrical transport in these systems be- 
cause one needs to contact tiny structures, a task that 
often turns out to be impossible. 

Among the techniques suitable for studying electronic 
transport in reduced dimensions, break junctions repre- 
sent an interesting option, because they provide an easy 
way to drive a current through only a few atoms. The 
measurement methods we used are based on the break- 
ing of a nanoscopic structure in a controlled manner while 
monitoring its resistance [10]. The setup is particularly 
stable since most of the structure is attached to the sub- 
strate and only the narrow bridge to be broken is sus- 
pended. It is then possible to mechanically stabilize the 
contact with a precision of a few picometres. We have 
slightly modified this procedure and decided not to sus- 
pend the bridge, because in ferromagnetic materials, any 
unsupported part is subjected to a magnetically induced 
distortion called magnetostriction. This effect results in 
a modification of the contact geometry when the mag- 
netization changes direction, which has been shown to 
significantly affect the contact resistance [11]. This dif- 
ference in the bridge fabrication is essential here as it 
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proves to be very efficient to reduce magneto-mechanical 
effects. Fig. 1 shows the result of a finite elements nu- 
merical simulation where a 20nm Fe layer is deposited on 
a kapton substrate with and without the under-etching 
procedure. Taking the physical parameters for the two 
materials and a saturation magnetostriction of polycrys- 
talline iron of —8ppm, the distortion of the Fe layer at 
the contact level goes from 15prn when the bridge is sus- 
pended to 1pm when the structure is attached to the 
kapton (no under-etching). We have checked that even 
in the tunnelling regime (where R changes exponentially 
with the gap), the effect is negligible. Direct forces due 
to stray fields are also found to be two orders of magni- 
tude smaller. In fact, any variation in the temperature of 
the electrodes would have a greater consequence than the 
magneto-mechanical effects. Therefore, the procedure in- 
sures that, in applying a saturating rotating field, the 
resistance change is only due to the intrinsic anisotropic 
"AMR" effect. 
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FIG. 1; Magnetostriction induced deformation of the nano- 
bridges in two geometries: (left) suspended and (right) at- 
tached to the kapton substrate (only half of the structures is 
represented - the other half is obtained by symmetry). The 
grey /colour scale represents the deformation which reaches 
15pm for the suspended bridge and 1pm when not under- 
etched (the field is along the bridge) . For visibility, the struc- 
ture with distortions magnified by 10000 is also shown. 

We have carried out an extensive set of measurements 
on Fe break junctions (similar results are also obtained 
for Ni and Co), where a 2.5 T magnetic field is rotated in 
the plane of the contact while the resistance is monitored. 
Interestingly, the resistance and the amplitude of its an- 
gular change did not significantly depend on the field 
magnitude between 0.5 and 2.5 T, where magnetostric- 
tion is expected to change by about 50%. This provides 
further experimental evidence that the effect can be ne- 
glected. In order to measure changes of resistance with 
magnetization direction, a high field was chosen to make 
sure the atomic contact is always in its saturated state 
(the demagnetization can be large in nanostructures). 

Figure 2 shows a representative set of curves in Fe at 
4.2 K. Interestingly, a behaviour qualitatively different 
from the cos^{0) dependence of the bulk, can be observed 
in the atomic contact regime. In the middle graph, it 
looks likely that one channel gets blocked when the field 
is along the contact, leading to a two-level conductance 
and an atomic-AMR (AAMR) effect of 21%. At slightly 
different values of conductance, both smooth sinusoidal 
variations as well as discrete jumps are observed (see top 
graph). This is to be expected when overlap changes are 
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FIG. 2: Variation of resistance as a 2.5 T field is rotated in 
the plane of the contact. The bottom graph is obtained in the 
first stages of pulling the bridge (as the nanostructure is not 
yet broken). There, the AMR is close to the bulk value, i.e. 
around 1%. In the atomic contact regime, at a conductance 
of 3e^//i in the middle graph the AMR behaviour is close to 
a two level effect reaching 3.5e^/h when the magnetization 
is perpendicular to the contact. In the top graph, where the 
conductance is close to 2e^ /h, the effect is intermediate be- 
tween the two behaviours and reaches 75% 



not sufficient to completely close a channel, but enough 
to change their transmittance. This general behaviour is 
consistent with what is known to happen for non mag- 
netic break junctions when orbital overlap is varied by 
mechanical deformation of the contacts [12]. There, the- 
oretical calculations have shown that both sharp jumps 
and smooth variations of the conductance can be ex- 
plained by considering the details of orbitals overlap [13]. 

Even more surprising, is the effect measured in tun- 
nelling and shown in figure 3. In this regime, charge 
carriers jump from one electrode to another through a 
very narrow vacuum gap (a few Aat most in break junc- 
tions). There, the evanescent wave functions still have 
a strong atomic orbital character from which they can 
inherit the spin-orbit coupling properties. Moreover, be- 
cause of the exponential decrease of the wave intensity 
with distance and the corresponding exponential increase 
of resistance, it is likely that any change in the shape of 
the orbitals could have a more significant effect. This 
is indeed what we experimentally observe. Interestingly, 
the two measurements shown in Fig. 3 correspond to elec- 
trons tunnelling through gaps of different sizes (as mea- 
sured in our setup) and similar resistance. Hence, the 
tunnelling cross section must differ which implies that 
electrons tunnel from sharp tips and small gap or flat 



110 
105 ■ 
100 ■ 









-30 30 60 90 120 150 180 210 240 

Angle (degree) 

FIG. 3; Measured resistance variation in Fe atomic contacts 
in the tunnelling regime when a 2.5 T field is rotated in-plane. 
The large 100% AMR efi'ect is obtained at a shorter gap value 
than that for the 20% effect of similar resistance. 
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tips and large gap. The atomic orbital character of the 
evanescent waves is obviously stronger in the first case, 
and one can expect the tunnelling-AMR (TAMR) to be 
larger, as measured. Moreover, because evanescent d or- 
bitals are generally shorter range than s orbitals their 
contribution should be larger for narrow gaps. 
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FIG. 4: Amplitude of the AMR effect in contacts of different 
resistances. The dotted line is a guide to the eye underlining 
the steep increase as the contact reaches atomic sizes. The 
vertical line is at /h and represents a rough separation be- 
tween tunnelling and atomic contact regimes. 

The measured ratio of high to low resistances (when 
rotating the field) for a Fe sample in different atomic 
configurations are gathered in figure 4. This plot is in- 
structive because it shows that the AMR effect increases 
steeply as the contact reaches atomic dimensions. For 
higher resistances (above e^/h), when in the tunnelling 
regime, we observe a significant scatter in the amplitude 
of the TAMR where values around 100% can be achieved 
but effects as low as 20% can also be found. As explained 
above, this can be understood because the atomic orbital 
character of the evanescent waves depends on the exact 



FIG. 5: Calculated band structure for a Fe mono-atomic chain 
(see top right inset) magnetized in the parallel b) and per- 
pendicular c) directions. In a) the band structure for a par- 
allel (full black circles) and perpendicular magnetization (red 
empty triangles) is shown at a larger scale around the Fermi 
level. 



atomic configuration, which can vary appreciably at con- 
stant resistance. 

Theoretically, a considerable amount of work has been 
devoted to metallic atomic contacts, but magnetism has 
seldom been considered. Most relevant works have stud- 
ied the resistance generated by a "magnetic domain wall" 
on the contact [3] and very recently, the contribution 
from the exchange splitting has been found to be rather 
small [4] , the dominating effect being instead the orbital 
nature of the conduction electrons [5] . In order to under- 
stand the origin of A AMR effect, we have performed ab- 
initio calculations to determine the changes of the elec- 
tronic band structures with the spins direction. In an ef- 
fort to give a pedagogical illustration of the phenomenon, 
we consider an ideal atomic chain of Fe atoms. The elec- 
tronic structure of the wire is obtained using the pseudo- 
potential plane-wave method implemented in the PWscf 
package [14] which, in its last version, allows to include 
spin-orbit interactions. Due to the use of plane waves, 
the system considered is in fact a periodic array of atomic 
chains, for which the distance between two wires is large 
enough (15 A) to avoid interactions. We have first carried 
out ultrasoft pseudo-potential calculations without spin- 
orbit in the Generalized Gradient Approximation (GGA) 
to obtain an equilibrium spacing a of the atoms in the 
wire of 2.27 Aand a magnetic moment of 3.3 /is per atom. 
Then using a fully relativistic (i.e solving the Dirac equa- 
tion for an atom) ultrasoft pseudo-potential including 
spin-orbit coupling [15] in the Local Density Approxima- 
tion (LDA) we have calculated the electronic structure 
of the wires with a lattice spacing of a. Fig. 5 shows the 
band structure obtained for the magnetization parallel 
and perpendicular to the chain. Significant changes near 
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the Fermi level are observed, which stem from the degen- 
eracies induced by spin-orbit coupling. The interesting 
bands are the weakly dispersive, so called Slater Koster 
ddS bands, which have a pure d^y (and dx2_y2) character 
since they couple neither to s and p nor to other d states 
with different symmetry. It is clear from Fig. 5 that these 
bands split by about 0.10 to 0.12 eV when the magne- 
tization is rotated from perpendicular to parallel to the 
atomic wire. It can actually be shown in a tight-binding 
model [16] that the splitting is equal to 2A, which leads 
to a value of A = 50 — 60meV. Because these states cross 
the Fermi level, they will play an important role in defin- 
ing conduction channels with some transmittance. Since 
they are also very sensitive to the magnetic orientation, 
a magnetoresistance effect can be expected. Indeed Fig. 
5 shows that when the magnetization is perpendicular to 
the chain the two degenerate ddS bands are crossing the 
Fermi level while in the parallel case their splitting 
pushes the upper band to higher energies where it hardly 
touches Ef. This results in a significant electronic trans- 
fer which almost empties the upper band. Therefore we 
can expect a better conductivity when the magnetization 
is perpendicular to the chain than when it is parallel, 
in good agreement with the experimental findings. Of 
course, a direct inference of the conductance from the 
band structure topology is not possible, but these quan- 
tities are intimately linked. What can then safely be said 
here is that the transmission of some conduction chan- 
nels will surely change with the magnetization direction. 
Beyond this simple atomic chain model, we have car- 
ried out preliminary tight-binding calculations in a more 
complex geometry of two pyramids in contact on their 
apex. The result is that the local density of states at the 
atoms forming the contact region is also modified when 
changing the orientation of the magnetization, although 
a bit less significantly than for the 1-D wire. Prom these 



examples, we note that the change in the energy level 
distribution is complex which could induce a positive or 
negative AAMR effect, rather sensitive to the atomic, 
magnetic, and electronic structure at the constriction re- 
gion. Hence, one can expect the AAMR to depend on the 
exact atomic geometry as well as the potential difference 
across the contact. One can even envision sign changes 
in the AMR with applied voltage. 

In conclusion, we have shown that rotating a saturat- 
ing field in an iron atomic contact leads to significant 
resistance changes. The effect is due to a spin-orbit cou- 
pling induced modification of orbitals overlap. This is 
supported by ab-initio calculations showing clear changes 
of the electronic band structure when the orientation of 
the magnetization in the samples is rotated by 90 °. The 
AMR effect also exists in the tunnelling regime where 
evanescent orbitals keep an atomic character. The mea- 
sured amplitude of the effect opens the possibility of us- 
ing it in field sensors or magnetic memories of atomic 
size. We find that the AMR effects (both AAMR and 
TAMR) are actually larger than the spin scattering ef- 
fects obtained when the magnetizations are opposite on 
the two sides of the atomic contact. 
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